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We study multiple helical nanowires in proximity to a common mesoscopic superconducting island, 
where Majorana fermion bound states are formed. We show that a weak finite charging energy of 
the center island may dramatically affect the low-energy behavior of the system. While for strong 
charging interactions, the junction decouples the connecting wires, interactions lower than a non- 
universal threshold may trigger the flow towards an exotic Kondo flxed point. In either case, the 
ideally Andreev reflecting flxed point characteristic for inflnite capacitance (grounded) devices gets 
destabilized by interactions. 



PACS numbers: 71.10.Pm, 73.23.-b, 74.50. +r 

Introduction. — Electronic transport through topologi- 
cal insulators |T] or superconductors |5] has come into 
the limelight of condensed-matter physics over the past 
few years. In particular, understanding the physics of 
localized Majorana bound states, generically expected 
near the ends of one-dimensional (ID) topological super- 
conductor wires |3jt5j, is crucial for exploiting the non- 
Abelian statistics of Majorana fermions in topological 
quantum computation schemes [6'-^. When a grounded 
Majorana nanowire — such as InSb and InAs wires sub- 
ject to a Zeeman field and proximity coupling to an s- 
wave superconductor |H1 [TU] — is weakly contacted by a 
normal metal, the presence of a Majorana state should 
reflect in a conductance peak |llH15j . and signatures of 
this type were indeed observed experimentally |16fll9| . 

In this paper, we discuss novel transport phenom- 
ena caused by Coulomb interactions in devices compris- 
ing Majorana wires contacted to leads. Our setup is 
sketched in Fig. [T] For N nanowires proximity-coupled 
to the same mesoscopic superconducting island, there are 
2N Majorana fermions, jj = 7j, with anticommutator 
{ijilk} = The island connects to the jth lead by 
tunneling through the Majorana state 7j with coupling 
strength tj ; all other coupling mechanisms are irrelevant 
[20J. Coulomb interactions affect the system in two dis- 
tinct ways: (i) The combination of repulsive interactions, 
spin-orbit coupling and Zeeman field is expected to turn 
each of the M < 2N connecting leads into a spinless 
(helical) Luttinger liquid (LL) characterized by an inter- 
action constant g < 1 |21fl23| . (ii) Weak intra-island in- 
teractions do not compromise the integrity of individual 
Majoranas |24fl26] . However, they introduce correlations 
between these states, and thereby correlations between 
the connecting wires. It stands to reason that the op- 
tion to generate inter-wire coupling on mesoscopic scales, 
i.e., independently of microscopic single-particle tunnel- 
ing between distinct terminals |27j i28j , will be crucial for 
quantum computation (and other) applications. 

Our main observations are summarized as follows. It 
turns out that Ec, no matter how weak, destabilizes the 



fully Andreev refiecting fixed point characterizing the 
grounded system |B] . Thinking of tunneling events from 
and into lead j in terms of particles and anti-particles, 
charging generates a strong confinement force operational 
at time scales r > E^^. At larger time scales, tun- 
neling events are bound into overall charge neutrality 
preserving 'dipoles', describing the near instantaneous 
(at scales E~^) transmission of charge from lead j to 
lead k ^ j. The formation of these objects terminates 
the up-renormalization of coupling constants prevalent 
in the 'asymptotic freedom' regime r < E^^. At scales 
T > E^^, the effective tunneling strengths are subject to 
a competition of downward renormalization due to fluctu- 
ations of individual dipoles, and upward renormalization 
due to dipole-dipole interactions. The balance of these 
two mechanisms deflnes a repulsive flxed point. A*, sep- 
arating a flow towards the decoupled dot (zero coupling) 
from a flow towards an exotic Kondo regime, generalizing 
the M = 3 topological Kondo effect discussed in Ref . [29] . 
Whether the system ends up in the trivial weak-coupling 
limit, or in the Kondo regime, depends on the bare cou- 
pling constants as much as on the value of Ec'. For weak 
charging, the initial growth of the bare charge coupling 
constants persists up to larger time scales, which gives 
the system an enhanced chance to reach couplings > A* 
before the conflnement regime takes over. We flnally note 
that in supporting a novel type of Kondo flow, the present 
system differs strongly from the M = 2 Majorana single- 
charge transistor |30ff52] as well as from other types of 
previously studied LL junctions |33H35| . 

Model. — We start by deriving the low-energy effec- 
tive theory for the setup in Fig. [l] The Hamiltonian, 
H = Hc + Ht + Hi, contains a piece He describing the is- 
land, the tunnel Hamiltonian Hf, and Hi for the LL leads 
(we put h = ks = 1 below). For each nanowire (a — 
1, . . . , N), we have two spatially well separated Majorana 
fermions, 720-1 and 720,. It is useful to deflne the nonlo- 
cal auxiliary fermion operators da — (720-1 + ij2a)/V2, 
where the total Majorana occupation number operator 
is h = ^a'^a^a- Assuming that the proximity gap ex- 
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Figure 1. (Color online) Schematic setup of the multi- 
terminal Coulomb-Majorana junction. A'' helical nanowires 
(here, A'^ = 3) connect to a floating mesoscopic superconduc- 
tor (shown as gray box) with charging energy Ec- Majorana 
fermions 7j (filled circles) exist near the ends of proximity- 
coupled wire segments. The remaining wire parts (away from 
the island) act as LL leads connected to the island by tunnel 
couplings tj. Majorana fermions representing Klein factors 
are shown as open circles. The case of M < 2N leads (here, 
M = 5) is included by putting one or several tj = 0. Direct 
tunneling between the 7j is neglected. 



ceeds all other energy scales, the island is fully described 
in terms of the and the Cooper pair number opera- 
tor, Nc, or conjugate phase, ip, where [ip,2Nc] = i. Since 
both, the Majoranas, and the superconducting phase are 
zero-energy modes, the island Hamiltonian is solely due 
to charging effects. 



(1) 



where rig is a dimensionless gate voltage parameter. The 
semi-infinite LL leads, with tunnel contacts at x = 0, are 
described by dual bosonic fields, 4>j{x) and 6j{x) |21| . 
where 



H, 
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27r ^ 



dx [gid^cb.y + g'\d,e^y 
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with Fermi velocity vp- The fermion annihilation op- 
erator for a right- or left-mover reads ipj.R/Lix) — 
^-i/2^^^i{<pj±ej) ^ where a is a short-distance cutoff and 
the r]j are additional Majorana fermions with {rjj, 7]^} — 
6jk [22]. These 'Klein factors' ensure fermion anticom- 
mutation relations between different leads and play an 
important role below. Since Hi describes decoupled leads 
with perfect normal reflection at a; = 0, we impose open 
boundary conditions, ipj^^iO) — 'ipj^f({0), pinning all fields 
Oj{0). A lead fermion operator near the respective tunnel 
contact is thus given by '^j — a~^^^r]je'^'^^^^\ Finally, the 
tunnel Hamiltonian reads [32] 

= E^^^*l (%/2] + (3) 
j 

where fj = (a/2)^''^ij. Tunneling either leaves no trace in 
the condensate (the first term) , or must involve the oper- 



ators e^^*''', raising or lowering the Cooper pair number 
Nc by one unit. 

Effective phase action. — Next we derive an equivalent 
Euclidean phase action, S, which will allow us to better 
understand the physics. We start by integrating over the 
harmonic LL bosons away from x ~ 0, retaining only 
$j(r) = (/)j(0,t). This generates the action piece Si ~ 
1§ n \^n\\^ji^n)\'^, with Matsubara frequencies Un — 
2TmT for temperature T (n G Z). Next, by virtue of a 
Hubbard-Stratonovich transformation to the condensate 
phase field (^(t), the charging term in Eq. ([T|) leads to 
the contribution 



Sr = 2'KiWn„ 
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The phase field has the periodicity (^(r -I- l/T) = ^{t) + 
2ttW with winding number g Z. However, below we 
mainly focus on the W = sector, which is appropriate 
for weak-to-intermediate charging energy. Noting that 
there also is a free fermion piece, Sj = / drdada + 
I J2j J drrjjrij, we now perform a gauge transformation, 
da{T) — i> e^*'^^'^^(ia(T), in order to eliminate the ipfi term 
in Eq. ([4|. It is then convenient to switch back from 
the d fermions to Majoranas, 72a-i = (^a + '^a)/'\/2 and 
72q — ~i{dc, — (iJ,)/\/2. After these steps, the tunnel 
action follows as 



St — E] tj / dT{—2irijjj) sin($j + ip). 



(5) 



At this stage it is useful to note that for each j, the 
Klein factor rjj and the Majorana 7^ can be fused to 
form the auxiliary fermion fj = {rjj — ijj)/^/2, where 
—2irjjjj — 2/J/j — 1 = (Tj = ±. The resulting fermion 

Hamiltonian conserves fjfj, and thus also each parity 
number aj. Hence Uj = ± is not a dynamical variable 
but just has to be summed over in the functional inte- 
gral The sole effect of this summation is to eliminate 
contributions of odd order in the tunneling amplitudes tj 
to the perturbative expansion of the functional integral, 
while even-order contributions remain unaffected. Keep- 
ing only even orders, we will drop the then immaterial 
presence of the cr^'s throughout. This Klein-Majorana 
fusion procedure implies an enormous technical simpli- 
fication, since we are now left with an effective phase 
action only. Before writing down this action, we shift 
$j — > ^j ~ if, which removes ip from the tunnel action 
([5]). Performing the remaining Gaussian functional inte- 
gral over ip, after unitary transformation to the discrete 
Fourier basis {q ^ 0, . . . , 2N — 1), 
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we arrive at the effective phase action 
Tg sr^ \ujn\ 



S = 



(7) 



with the energy scale e = AgNEc/ir. We stress that the 
charging energy only affects the zero mode $o, which 
effectively becomes free at low frequencies, |a;„| ^ e. 

Particle analogy and renormalization. — We next turn 
to discussing the perturbative expansion in the tunnel 
couplings tj. A very instructive analogy follows by in- 
terpreting 0^{t) = e^**^^'^-' as particles ('quarks') and 
antiparticles living on the imaginary-time axis. Our 
particles carry a 'flavor' index (j = f,...,2iV) and an 
(imaginary- valued) effective charge —itj/2. We study the 
properties of the ensuing interacting particle gas by stan- 
dard renormalization group (RG) methods [53^, "371 . In an 
RG step, we integrate over fast modes with frequen- 
cies in the shell A/6 < |aj„| < A, with rescaling parame- 
ter & > 1 and a high-energy cutoff A of the order of the 
proximity gap. The RG step is completed by rescaling 
all frequencies, uj buj, such that A stays invariant. For 
small tj, the particle density is low and different charges 
(i.e., the tj) renormalize independently. We first RG- 
integrate out modes in the shell e < |a;„| < A, where 
the zero mode $0 is not significantly affected by Ec, see 
Eq. ([7|. Some algebra yields the net scaling dimension 
1 — l/(2g) for tj, which are therefore RG-relevant cou- 
plings for g > 1/2. This signals a flow towards the res- 
onant Andreev reflection fixed point [20_. In the lead- 
non-interacting case, g = I, the scaling dimension 1/2 
represents the naive dimension of a fermion-Majorana 
scattering operator. After the RG integration over the 
shell e < |a;| < A, the renormalized couplings are then 
given by ^ = tj (A/e)i-i/(29). 

Proceeding to lower frequency scales, |a;| < e, the 
charging energy renders the zero-mode action non- 
dissipative, S[^q] ~ JJ^ |l>o(w„)|2, see Eq. ^. 

The integration over $0 generates a strong 'confinement' 
potential linear in the time separation between particle 
creation, Oj'{T) ^ q+'^'^o{'^)/V2N ^ ^^^^ particle annihila- 
tion events, {t') ^ e-**o(-r )/V2N q£ arbitrary fiavor 
index k ^ j. (For k = j, particle-anti-particle annihi- 
lation occurs - the inconsequential virtual tunneling of 
a particle to and fro the same lead, cf. Fig. [2] upper in- 
set.) At time resolutions lower than E^^, the relevant 
excitations of the theory then are quark-antiquark pairs 
('mesons') [55] . 



(0+(r)0-(r'))o 



"^'1 0,fc(r), 
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RG equations in dipole regime. — We now consider the 
effective low-energy theory at resolutions EcT > 1, where 



we are dealing with a gas of dipoles, Ojk, with couplings 
^jk, j 7^ k. The 'bare value' of these couplings is given 
by Xfj « tf k^^'>/Ec > 0, where the factor E^^ is due to 
a time integration over the distance between particle and 
antiparticle. Physically, these couplings bundle the effect 
of in-tunneling from lead j into a virtual on-dot state of 
longevity ~ E^^-, followed by out-tunneling into lead k. 
We then continue the frequency-shell RG integration for 
mesons, going up to one-loop accuracy. Some algebra 
yields the RG equations 



d\nb 



= (1-5-') A,* 



K 

e: 



(9) 



where k is a non-universal constant of order unity. The 
first term describes the standard power-law suppression 
of the tunneling density of states in/out of the LL leads 
|21| and implies a suppression of the Xjk- This term 
is fought by the positive second contribution, which de- 
scribes the effect of dipole fusion [(j, to) -|- (to, k) — >■ (j, k)] 
due to dipole-dipole interactions (cf. Fig. [2] lower inset.) 
Equation (|9| suggests the existence of an RG-unstable 
fixed point with isotropic couplings. 



Xjk — A — 



- 1 



k(M-2) 



Er. 
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Discussion. — When decreasing the effective frequency 
scale a;, typical couplings, t and A, show the RG flow il- 
lustrated in Fig. [2] During the initial RG flow (down 
to oj w Ec), particles are asymptotically free, imply- 
ing a power-law increase of a typical tunnel coupling, 
t oc Lu^n , with the time-like variable w . We now 
compare the value A^^^ = {t^^^)'^/Ec Ec^^^^^, reached 
at the end of the asymptotic freedom phase (w w Ec), 
to the fixed-point value A* ^ Ec in Eq. (10 1. For suffi- 
ciently large Ec [H], the fixed point cannot be reached, 
and therefore the w — >■ stable fixed point describes 
a completely decoupled LL junction. For sufficiently 
small Ec, however, the fluctuations generated during the 
asymptotic freedom phase tip the balance, A*^^-* > A*. 
A straightforward expansion near the fixed point then 
shows that the isotropic baseline A > A* of the cou- 
plings fiows to large values with dimension g^^ — 1 > 0, 
while deviations between the couplings are RG irrele- 
vant, i.e., the flow is towards an isotropic conflgura- 
tion Xjk — X. At a characteristic 'Kondo temperature', 

Tk ~ Ec exp ^(1) i^fj^i _2) ^ ■ the coupling A begins to di- 
verge, and our perturbative expansion is no longer appli- 
cable. A comprehensive discussion of the ensuing strong- 
coupling regime is beyond the scope of the present paper. 
However, by analogy to the M = 3 topological Kondo 
fixed point [29], we expect that A — )■ cx) is the only fixed 
point besides A = 0, A*. 

Conductance matrix. — The two-stage scenario outlined 
above, with its branching into either a strong-coupling or 
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reveals that the physics of the system is determined by 
asymptotic freedom at short time scales (r < E~'^) and 
confinement at long time scales (t > E~^) correspond- 
ing, respectively, to independent dot-lead tunneling, and 
virtual co-tunneling. For sufficiently weak Ec, we find 
that a strong-coupling Kondo fixed point is approached, 
while otherwise the junction describes decoupled leads at 
T — Q. These predictions can be tested using the tem- 
perature dependence of the conductance. 

We thank P. Sodano and A. Zazunov for discussions. 
This work was supported by the SFB TR 12 of the DFG. 
Note added: During the preparation of this manuscript, 
we learned of independent related work by Beri |40| . 
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Figure 2. (Color online) Schematic RG flow of typical cou- 
plings, t and A, vs the time-like RG parameter oi"^, and il- 
lustration of the system excitations. At short times, indi- 
vidual in- and out-tunneling events into difi'erent leads domi- 
nate (indicated by dots and circles of difi^erent color). These 
'quarks' are asymptotically free, and their fluctuations gener- 
ate a power-law increase of t. For w < Ea, confinement leads 
to bound quark-antiquark dipoles (mesons) for j / fc, or an- 
nihilation for j = k. The dipole coupling. A, experiences a 
downward renormalization due to individual fluctuations and 
an upward renormalization due to dipole-dipole fusion events 
(see also right lower inset). For sufBciently large Ec — Ed, A 
does not reach the fixed-point value ( 10 1, the first mechanism 



dominates, and hence A flows to the decoupled fixed point as 
oj — > 0. For small Ec = Ec2, however, the increase of t during 
the asymptotic freedom phase brings A beyond the balance 
point A* (see text), and the RG fiow proceeds towards the 
strong-coupling topological Kondo fixed point. 



a decoupled low-frequency limit, will bear consequences 
for all physical observables. We here briefly discuss the 
resulting temperature dependence of the linear conduc- 
tance Gjk between different leads j and k. For high 
T > Ec, asymptotic freedom causes the power-law scal- 
ing Gjk ^ r^2+i/£(^ while for T — > 0, we either have a 
vanishing conductance, Gjk ^ T^^+^/f , if the decoupled 
fixed point is approached, or the Kondo-type flow ap- 
proaching the unitary limit of an isotropically hybridized 
junction, Gjk = ^jj- The latter case constitutes the 
M-terminal generalization of the teleportation scenario 
discussed for M = 2 by Fu [31] . It is worth stressing that 
for both stable fixed points (Kondo or decoupled), reso- 
nant Andreev refiection is unstable, and hence arbitrary 
Ec destroy the corresponding fixed point. 

Conclusions. — In this paper, we have formulated and 
studied the problem of junctions of Majorana wires meet- 
ing on a superconducting island with charging energy Ec- 
We also included correlations in the leads, since these 
typically are ID nanowires themselves. Our RG analysis 
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